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We compute the nonzero temperature free energy up to the order g 6 ln(l /g) in the coupling constant 
for vector like SU(N) gauge theories featuring matter transforming according to different representa- 
tions of the underlying gauge group. The number of matter fields, i.e. flavors, is arranged in such 
a way that the theory develops a perturbative stable infrared fixed point at zero temperature. Due 
to large distance conformality we trade the coupling constant with its fixed point value and define a 
reduced free energy which depends only on the number of flavors, colors and matter representation. 

We show that the reduced free energy changes sign, at the second, fifth and sixth order in the 
coupling, when decreasing the number of flavors from the upper end of the conformal window. If the 
change in sign is interpreted as signal of an instability of the system then we infer a critical number 
of flavors. Surprisingly this number, if computed to the order g 2 , agrees with previous predictions 
for the lower boundary of the conformal window for nonsupersymmetric gauge theories. The higher 
order results tend to predict a higher number of critical flavors. These are universal properties, i.e. 
they are independent on the specific matter representation. 
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I. INTRODUCTION 

Non-Abelian gauge theories are expected to exist in 
a number of different phases which can be classified 
according to the force measured between two static 
sources. The knowledge of this phase diagram is rel- 
evant for the construction of extensions of the Standard 
Model (SM) that invoke dynamical electroweak symme- 
try breaking 1 1 , 2 1 . An up-to-date review is ]3) while 
earlier reviews are |4, 5 J. The phase diagram is also use- 
ful in providing ultraviolet completions of unparticle 1 6 1 
models [7 . 8 [ and it has been investigated recently using 
different analytical methods I19WT91. 

Here we wish to understand, in a rigorous way, the 
dynamics of gauge theories lying in the conformal win- 
dow at nonzero temperature. The physical applications 
are numerous ranging from the above mentioned mod- 
els of dynamical electroweak symmetry breaking to cos- 
mology [20-24]. Thermodynamical properties of these 
gauge theories were also investigated in the literature 
using holographic models [17ll25ll2"6] . 

Our starting point are asymptotically free vector-like 
gauge theories near the Banks-Zaks infrared stable fixed 
point 1 27]. The presence of such a perturbative fixed 
point allows a controllable computation of the free en- 
ergy for these theories which we carry till order g 6 in the 
gauge coupling. In fact we expect perturbation theory 
to work, in principle, for any temperature in contrast 
to the case of a confining theory for which perturbation 
theory is limited to asymptotically large temperatures 
compared to its intrinsic scale. The absence of an intrin- 
sic scale in the theory is evident in having a free energy 
directly proportional to the fourth power of the tempera- 
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ture, for any temperature. We will trade the value of the 
coupling with its value at the infrared fixed point turn- 
ing the coefficient of T 4 into an algebraic expression of 
the number of flavors, colors and matter representation, 
encoding a great deal of information of the underlying 
gauge theory. 

Despite the fact that we can keep the coupling con- 
stant small it is a fact that perturbation theory breaks 
down, at finite temperature, due to the loss of analytic- 
ity in the coupling associated to the presence of infrared 
singularities. For the free energy this problem sets in at 
0(g 6 ), or four-loop order|28 ,29\. At best one can assume 
that the free energy is computable to 0(g 6 ), though not 
via loop diagrams. The highest order one can achieve 
using Feynman diagrams is O (g b ln(l / g)J and was re- 
cently determined in [30 1 . We adapt their results for the 
case of gauge theories featuring large distance confor- 
mality while generalizing the discussion to any matter 
representation, different number of colors and flavors. 

We discover a number of surprising features when 
plotting the free energy, for a given matter representa- 
tion, as function of the number of flavors: i) There is 
a change in sign of the free energy at a critical number 
of flavors whose value depends on the representation it 
belongs and the order to which the computations were 
carried, ii)This number is smaller than the one for which 
asymptotic freedom is lost. 

It is tempting to identify it with the critical number of 
flavors below which, at zero temperature, conformality 
is lost. The obvious caveat is that as we decrease the 
number of flavors away from the point when asymp- 
totic freedom is lost perturbation theory ceases to be reli- 
able and therefore we interpret this phenomenon only as 
strong indication that the finite temperature free energy 
is aware of the nontrivial underlying gauge dynamics. 

Another amusing feature is that at the two-loops level, 
when the results are scheme independent, the change 
in sign of the free energy occurs for a given number 
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of flavors which is surprisingly close to the one pre- 
dicted using the Schwinger-Dyson results IflOl as well 
as the Ryttov-Sannino /3 function El . This value be- 
comes larger when going to the fifth and sixth order in 
the coupling. 



II. REVIEW OF THE FREE ENERGY COMPUTATION 

The perturbative free energy at finite temperatures 
was computed to order O (g 6 ln(l/g)) in p0| and in the 

MS scheme. The result was derived using an effective 
field theory approach utilizing matching of coefficients 
in the effective theory expression with the dimension- 
ally regularized perturbative expansion. An in depth 
presentation of the method can be found in llSTl , where 
the order 0(g 5 ) was recomputed. In order to present 
the results in a relatively self-contained way we briefly 
review the method here. 

The perturbative free energy for a massless asymp- 
totically free gauge theory receives contributions from 
the following three mass scales: 2nT, gT and g 2 T. They 
are respectively the particle momentum in the plasma, 
the onset of the color-electric (Debye) screening and fi- 
nally the onset of color-magnetic screening. The idea is 
then to construct an effective field theory that reproduce 
static observables at the different scales. This is done by 
the method of dimensional reduction, where the static 
properties of a 3 + 1 -dimensional field theory at high 
temperatures are expressed in terms of an effective field 
theory in 3 space dimensions 129 l32l . 

The free energy density expressing the static equilib- 
rium properties of the plasma is given by the usual log- 
arithm of the partition function: 



T 

F = --lnZ, 



z 



j DApDxpDipexpi- d% j d d xA, 



(1) 

(2) 
(3) 



As it is customary to introduce finite temperature by 
euclideanizing the time dimension t = -h and compact- 
ifying it with a period |3 = 1/T, with T the temperature, 
and the bosonic (fermionic) fields respecting periodic 
(antiperiodic) boundary conditions. In regimes when 
the gauge coupling, g is small the free energy density 
can be computed perturbatively. 

On the other hand, by dimensional reduction one can 
compute the free energy density using an effective field 
theory in 3 space dimensions, by seperating the electro- 
and magnetostatic parts of the Lagrangian. The free 
energy density is expressed as 



F - T [f E (T, g; A E ) + /mK, gE, A®, . . . ; A E , A M ) 
+/g(£m/..-;Am)], 



(4) 



where the effective free energy densities /e//m and fc 
represents the contributions from the three scales, i.e. 
/e gives the contribution from the momentum scale by 
effectively integrating out the fermions and the high mo- 
mentum degrees of freedom down to the scale Ae corre- 
sponding to a distance of order l/(gT). At greater dis- 
tances the fields are replaced by electrostatic and magne- 
tostatic gauge fields Aq(x) and A"(x), which are propor- 
tional to the zero-frequency modes of the gauge fields 
A^(t, x) that are the only fields able to propagate over 
such distances 1 33 1 . /e is then the normalization for this 
transition, i.e. 



Z = e 



-/e(A e )V 



DA"(x)exp 



(-1 



d d x£ E ) 



(5) 



where the Lagrangian is now an effective electrostatic 
Lagrangian: 



-TrFif + Jr[Di,A ] z + m 2 E TrA 2 



+ a£> [Tr(A 2 )f + AfTcA* + 5£ E . (6) 



<5-£e represents higher order interaction terms which con- 
tributes beyond the order 0(g 6 ). Shorthand notation has 
been used for the gauge fields, i.e A f , = T"A", and the 
corresponding gauge coupling is denoted by gE- This 
Lagrangian defines /m which is dependent on the still 
lower momentum scale Am, corresponding to the dis- 
tance l/(g 2 T). Again at greater distances, only the mag- 
netostatic gauge fields play a role: 



Z = e 



-/e(Ae)V„-/m(A e ,Am)V 



f 



£>A?(x)exp - d d x£ 



(-;• 



-M 



with the effective magnetostatic Lagrangian 



£m = -Tr(F?) + 5£ M/ 



(7) 



containing the gauge coupling g M . Note that the di- 
rect identification of the normalization functions /e and 
/m with the free energy densities of the respective La- 
grangians is strictly true only when using dimensional 
regularization to cut off the ultraviolet and infrared di- 
vergences in the perturbation expansions. In the same 
sense, fc is identified with the integral expression in the 
above partition function. 

The Lagrangian |7]| defines a confining theory when 
the higher order terms are not considered and is thus 
non-perturbative. However, fc can still be expressed as 
a power series in g [33] . The leading order is proportional 
to (g 2 T) 3 and the coefficient can be determined by lattice 
computations only. However, the logarithmic ultravio- 
let divergence coming from the scale Am can be evalu- 
ated exactly by matching the coefficient with that in /m, 
since the expression must be scale invariant. In this way 
one proceeds backwards and matches all coefficient with 
appropriate tuning to get rid of the somewhat arbitrary 
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momentum scales. Finally one ends up with an exact ex- 
pression for the free energy to order C>(£ 6 ln(l/£)), while 
the order 0{g 6 ) coefficient remains unknown and un- 
computable from perturbation methods. One also shifts 
the scale of the coupling constant from the dimensional 
regularization scale to an arbitrary renormalization scale 
H by using the renormalization group equation for the 
running of the coupling constant. 

We here express the leading order magnitudes only 
and refer to 1301 for the full result. 



^~g% A?~g% <? M ~?T. 



h 



Note that the perturbation expansion parameter in Xe 
is gi/rriE which is of order g. Thus the perturbation 
expansion by this method is an expansion in g rather 
than g 2 . 



III. HOT CONFORMAL FREE ENERGY @ 0(g 2 ) 

Our starting point is a generic asymptotically free 
gauge theory with Nf Dirac flavors transforming ac- 
cording to the representation r of the underlying gauge 
group. We will consider, to this order, also the case of 
N = 1 supersymmetric gauge theories for reasons which 
will become clearer shortly. 

To be more specific, throughout this paper we will 
consider matter transforming according to four different 
but single representations, i.e. the adjoint representa- 
tion (denoted G) under the gauge group SU{2), the fun- 
damental representation under the gauge groups SU(3) 
and SU{2), the two-index symmetric representation un- 
der the gauge group SU(3) and the two-index antisym- 
metric representation under the gauge group SU{A). 

The relevant group normalization factors are: 



TxlT r T h r ] = T[r]8 



ab 



TX = C 2 [r]l, 



(8) 



where T" is the fl-th group generator in the representation 
r and a = 1, ... , d[G]. We denote with d[r] the dimension 



of the representation. T[r] and C 2 [r 



are related via the 



identity C 2 [r]d[r] = T[r]d[G]. In table ij for completeness, 
we list the normalization used for the group factors in 
the different representations. We list in the last column 
also the number of colors which will be considered. The 
normalizations were taken from [12J. 
The f> function up to four-loop order 



_ft_ 
{An) 2 



ft 
{An) 4 



ft 

{Anf 



ft 

(47T) 1 



g 9 + 0{g n ), 
(9) 



was computed in [34 J. As for the free energy expres- 
sion the four-loop |8 function is also computed in the MS 
scheme, thus no ambiguities in the scheme-dependence 



TABLE I: Normalization of the relevant group factors 
for the representations used throughout this paper. 
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of the expressions arise. Only jSo and |3i are scheme- 
independent and read: 



ft = yC 2 [G] - \T[r]N f , (10) 

34 / 20 \ 

ft = y C 2 2 [G] - (y C 2 [G] + 4C 2 [r]j T[r]N f . (11) 



Asymptotic freedom is lost when the lowest order coef- 
ficient, jSo changes sign. This occurs at 



AF 11 C 2 [G] 
/ ■ 4 T[r] ' 



(12) 



For a given fermion representation, the second coeffi- 
cient, jSi is negative below this critical number of fla- 
vors and an infrared-stable fixed point develop which is 
known as the Banks-Zaks fixed point[27J. Such a theory 
display large distance conformality. 

The Banks-Zanks fixed point disappears when jSi 
changes sign. This occurs at: 



\7C 2 {G] C 2 [G] 
/ ■ 10C 2 [G] + 6C 2 [r] T[r] ' 



(13) 



We are now equipped to investigate the conformal free 
energy by starting with the nonsupersymmetric case to 
the order g 2 : 



d[G] [1 7 d[r] 



n 2 T 4 



N 



Of) 

{An) 2 



20 d[G] 

For the supersymmetric case we have, 
)Sg USY = 3C 2 [G] - 2T[r]N / , 
jSf SY = 6C^[G] - 4(C 2 [G] + 2C 2 [r])T[r]N f , 
leading to: 

3 C 2 [G] 



N AF 

^f,SUSY 



N 



Lost _ 



2 T[r] ' 
3C 2 [G] 



C 2 [G] 



f,susv 2C 2 [G]+4C 2 [r] T[r] ' 



(14) 

(15) 
(16) 

(17) 
(18) 
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We obtain for the supersymmetric free energy[35|: 
Fsusy d[G] 



7I 2 T 4 



24 L 



1 +2 dW] Nf 



6(c 2 [G]+6T[r]N f ) 



2 (f) 



(4tt) 2 



(19) 



To determine the free energy dependence on the num- 
ber of flavors and colors in the perturbative regime of 
the conformal window we replace the coupling constant 
with the Banks-Zaks fixed point value g* at two-loop 
order, given in appendix [A| The free energies read: 



d[G] 



T 2 T 4 



5 20 d[G] 



f 



(C 2 [G] + lT[r]N f ) (11C 2 [G] - 4T[r]N f ) 



34C 2 2 [G] - (20C 2 [G] 



12C 2 [r])T[r]N f 



SUSY 

tz 2 T 4 



d[G] 
24 



1 + 9 d[r] N 
1+2 d\G] f 



3 (C 2 [G] + 6T[r]N f ) (3C 2 [G] - 2T[r]N/) 



3C 2 [G] - 2(C 2 [G] + 2C 2 [r])T[r]N f 



We observe immediately that due to the conformal large 
distance nature of our theories the free energy depen- 
dence on the energy scale is only via the temperature 
which factors out leaving behind, as expected, a numer- 
ical factor containing information on the specific theory 
studied. These coefficients are universal, i.e. indepen- 
dent on renormalization schemes. 

We can now plot the free energies as function of num- 
ber of flavors for different number of colors and gauge 
theories. We choose to normalize our results to the free 
energy obtained, at order g 2 , when replacing the number 
of flavors with the one for which asymptotic freedom is 
lost for any given underlying gauge theory (defined as 
Fj in the plots). The results are shown in figure [T] A 
generic feature emerging from the plots is that there is 
always a critical number of flavors N f for which the nor- 
malized free energy vanishes. We interpret the change 
in the sign of the free energy as an indication of an insta- 
bility of the system which identify with the point where 
large distance conformality is lost. In table III] we com- 
pare this value with the expected critical number of fla- 
vors obtained using the Ryttov-Sannino |3 function N^ s 

(obtained for the parameter y = 1) as well as the one ob- 
tained via the ladder approximation and indicated with 

\jLadder 
f " 

The agreement among these numbers is surprisingly 
good as it can been from the column indicated by dis- 
crepancy defined as \Nf - N^ s \/N^ s . 

We still do not have a deep understanding of why 
these different methods agree so well among each other 
however we speculate that this agreement might be due 
to the fact that all these approaches make use of the two- 
loop universal coefficients of the jS function. 



TABLE II: Comparison of the different critical 
number of flavors obtained via the g 2 free energy 
(Nf), the Ryttov-Sannino f> function N^ s (obtained for 
the parameter y = 1), and the ladder approximation 

\jLadder 

7 ' 
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2.2 


0.00699 
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7.35712 


7.33 . . . 


0.00324 


8.104 
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2 


1.82739 


1.833... 


0.00324 


2.075 



Once we noticed such an agreement we asked our- 
selves: How about supersymmetry? We find, also in the su- 
persymmetric case, the existence of a critical number of 
flavors below which the free energy changes sign. In ta- 
blefT 



III 



[ we compare Nf with the critical number of flavors 
obtained using the supersymmetric all-orders jS function 
when setting to zero its numerator for both y — — 1 and 
y = — |. In this case we find a reasonable agreement 
when taking as critical number of flavors the one for 
which the anomalous dimension of the chiral superfield 
y is around -1/5. This is not the preferred value ob- 
tained from Seiberg's results[36| which, however, were 
tested via dualities only for the case of the fundamental 
representation. 



IV. CONFORMAL FREE ENERGY TO THE LAST 
PERTURBATIVE ORDER 

Having at hand a perturbative expansion it is natural 
to go beyond the g 2 order. To determine the free energy 
at any given order, in perturbation theory, we have con- 
sistently solved for the value of the coupling evaluated 
at the infrared fixed value and in the same renormaliza- 



TABLE III: Comparison of Nf for supersymmetric 
gauge theories with the critical number of flavors 
obtained using the supersymmetric all-orders f> 
function when setting to zero its numerator for both 
y — — 1 and y = — \ 
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N f 
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0.59 (0.042) 
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4.9794 


3.0 (5) 




0.77 (0.004) 
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1.2071 


0.75 (1.25) 




0.61 (0.034) 
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(a) S(J(3) with fundamental 
fermions. 



(b) SU(2) with fundamental 
fermions. 



(c) SU(3) with two-index 
symmetric fermions. 



F 
F 




F 
F 




(d) SU(4) with 2-index (e) SU(2) with adjoint fermions. 

antisymmetric fermions. 

FIG. 1: Normalized free energy in the conformal window for different theories. 



tion scheme. The expressions of the fixed point value 
of the coupling to the highest order computed here are 
given in the Appendix [A| 

The four-loop |3 function was computed in [34 J up to a 
normalization constant for the fourth order Casimir. In 
11341 the explicit expressions for all the coefficients were 
given for the fundamental and adjoint representations 
while we derive in Appendix [B] the expressions for any 
totally (anti)symmetric representation for SU(N), SO(N) 
and Sp(N) gauge groups. 

Beyond the g 3 order one notices the emergence of log- 
arithms of the ratio of the renormalization to the temper- 
ature scale. Since we assumed, in our computations, the 
temperature scale to be such that the gauge theory cou- 
pling constant, at zero temperature, has (quasi) reached 
the fixed point value it is therefore natural to evaluate 
the coupling at the renormalization scale point 2nT. We 
have, however, checked by direct evaluation of the free 
energy at the renormalization scale point of g 2 T that, due 
to the logarithmic dependence, the results are rather in- 
sensitive to the choice of the reference scale as it is clear 
from Fig. [2] where we show the results for fermions in 
the fundamental representation for the two choices of 
the renormalization scale point, 2nT (left-panel) and g 2 T 
(right-panel). 

It is for this reason that we show in Fig.|3]the results for 



the remaining theories evaluated at the scale 2nT. We 
observe the following universal behaviors: 

• The fee energy to the lowest interesting scheme- 
independent order in perturbation theory (i.e. g 2 ) 
changes sign at a critical number of flavors (Nf), 

• This critical value increases at the order g 5 and 
increases further at the order g 6 , 

• The free energy does not change sign if truncated 
at the order g 3 or g A . 



V. CONCLUSIONS 

We unveiled the finite temperature structure of gauge 
theories of fundamental interactions featuring a per- 
turbative infrared stable fixed point to the last com- 
putable order in perturbation theory. Differently from 
gauge theories assumed to generate a nonperturbative 
renormalization-invariant scale at zero temperature, like 
QCD, our results are perturbative in the entire energy 
range (i.e. for any choice of the temperature) since we can 
use as control parameter the number of flavors to tune 
the theory near the perturbative stable infrared fixed 
point. 
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We discovered a number of universal properties, i.e. 
independent on the matter representation and the su- 
persymmetric structure of the underlying gauge theory 
suggesting that asymptotically free gauge theories fea- 
turing large distance conformality share very similar dy- 
namics. 

If we were to take the point of view 1 30 1 that having ex- 
hausted the perturbative results we determined the full 
result for the free energy at nonzero temperature, we 
would then have discovered that there is a critical num- 
ber of flavors below which the free energy changes sign 
signaling the onset of an instability which we interpret 
as the end of the conformal window. 
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FIG. 2: Normalized free energy for different orders in g computed at the renormalizations scale [i with fermions 
in the fundamental representation of SU(3). 




(a) SU(2) with adjoint fermions. (b) SU(3) with two-index symmetric (c) SU(4) with two-index antisymmetric 

fermions. fermions. 

FIG. 3: Normalized free energy for different orders in g computed at the renormalizations scale /i = 2nT with 
fermions in different representation. The color-code is the same as in figure [2] 
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Appendix A: Banks-Zaks fixed points up to four-loops 

Here we give the exact expression for the Banks-Zaks 
infrared fixed point \27\ to different orders in g. 
The two-loop expression is: 



The fourth-order Casimir invariant is related to the 
symmetrized fourth-order trace: 



of 

471 



(Al) 



The three-loop expression is: 



a 
An 



(A2) 



The four-loop expression is: 

a* ^(l + zV3)-&i(l + jV3) + 2jS 2 & 2 



4k 
where 



h = I2/S1/S3 - 4$, 



12^ 3 &2 



(A3) 



108j8 /3| - : 



&2 = (36jSi/3 2 /3; 

+12^/3^3-^ 
with the |3,'s given in 



18j3ij3 2 j3 3 )3o + 27 J 8 2 lJ S 2 , + 4/S j6|jS 3 ) 3 , 



Appendix B: Generalization for the 4th order Casimir 

The general fourth-order Casimir invariants for sim- 
ple Lie groups were derived in ||37l |38] . We will use 
the results therein to generalize the expression for the 
four-loop jS function in [34] to any representation of the 
groups SU(N),SO(N),Sp(N). We must however keep 
track of the different normalization of the Killing form 
in the literature, thus we define an overall normalization 
constant, b as 



1xf a J bcd = bhb ab = rjI 2 [G]S 



nb 



(Bl) 



where f abc are the structure constants and h is the dual 
Coxeter number. The equation defines the second-order 
Casimir invariant, h with eigenvalue h[A as given in 
[37 \ and r\ = b/b' relates the results therein to the arbi- 
trary normalization b, with V chosen in [37.. 38 1 to be 
V = {2, 1, 2} for the groups (S!i(N), SO(N), Sp(N)}. In this 
paper, b is set to 1. In this appendix we are following 
the notation introduced by Okubo 1 37 1 when naming the 
second-order Casimir invariants, i.e. li[r\, correspond- 
ing to C2I/] in the more recent literature and the one used 
in the main text. 



-1 

df cd = - Tr |r fl T b T<T d ] , 



(B2) 



where the sum is over all permutations of the generator 
indices. We let T" be any representation of the generators 
for a simple Lie group. It is well-known that the fourth- 
order Casimir invariant is not unique, in fact the square 
of I2 is also a fourth-order Casimir invariant, leading to 
the amibiguity: 



I' 4 = h + C(I 2 ) 2 , 



(B3) 



where I' 4 defines a new fourth-order Casimir invariant, 
with C being an arbitrary constant. To get rid of the 
ambiguity one defines a modified Casimir invariant J$ 
with a specific metric found by Okubo and for short 
indicated with 5 abcd : 



h = r]^ abcd T a T b rT d . 



(B4) 



Then requiring the identity for irreducible representa- 
tions: 



Jabcd 



Tr |Y fl x b r : 'r rf J = r] 2 rf[r]/ 4 [?- 



(B5) 



where ]a{t[ is the eigenvalue of J4 in the representation r, 
it follows that J4 satisfies similar sum-rules as I2 and I3, 
hence it is the appropriate fourth-order Casimir invari- 
ant to work with ||37| . 

It now follows that for a general representation, one 
can write: 

df cd = c Y b abcd + X -C2{b ab b cd + b ac b bd + b ad b bc ), (B6) 

where c, are some constants dependent on the represen- 
tation, r. The two terms are orthogonal. We will only 
be concerned with the first term, as the second term was 
given for any representation in [ 34 1 . From |37] one finds 
that 



d[r]k[r] 



d[A]/ 4 [A](2 + d[G))' 



(B7) 



where A is the defining representation, which we will 
take to be the fundamental one. Then, contracting eq. 
dB6l with b abcd one finds: 



Oabcd 



4A]/ 4 [A](2 + d[G]). (B8) 



Hence, we derive that 

d[r\h[r} 



Jlfjibcdf. 



d[A]/ 4 [A] ' 

d[r]hM , 2 
4A]/ 4 [A] 



d[A]hW 
(2 + d[G]) 

d[A]hW 
b'\2 + d[G])' 



(B9) 
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Writing the expression in this form, we exactly get the 
definitions of the normalization constant I^r] and the 
traceless tensor d abcd used in the four-loop |3 function 
paper |34J (note that the normalization constant was de- 
fined without the tilde, but is used here in order not to 
confuse it with I4 [r] defined in |37|), i.e. 



SO(N): 



d abcd d, 



j r ! d[r]hlr] , 2 
d[A]JM] 



abed 



fe' 4 (2 + d[G])' 



(BIO) 
(Bll) 



As noted in |34|, d ahcd is representation-independent, and 
the contracted product can be written as 



d abcd d abcd [SU(N)] = 



d abcd d abai [SO(N)} = 



d abcd d abcd [Sp(N)] 



d[G](d[G] - 3)(d[G] - 8) 
16 • 6(2 + d[G]) ' 

d[G](d[G] - l)(d[G] - 3) 
12(2 + d[G]) ' 

d[G](d[G] - l)(d[G] - 3) 
16 ■ 12(2 + d[G]) 



(B12) 
(B13) 
(B14) 



where the fundamental representation is taken as the 
defining representation A, and with 





SU(N) 


SO(N) 


Sp(N) 


d[G] 


N 2 -l 


N(N - l)/2 


N(N + l)/2 



Correspondingly Ii[r] can be derived from |37|, where 
all invariants are given. For completion, we give the 
expressions here. 

Denote the fundamental representations with {Aj} cor- 
responding to completely antisymmetric tensor repre- 
sentations, while {AAi} are the completely symmetric 
representations in the sense of Young's tableaux, i.e: 



A, 



m 



AAi ~ EE ■ ■ • E 



To simplify the expressions, we define: 

Cj=N(N + l)-6j(N-f), 
K k = N(N - 1) + 6k(N + k) 

Then, we find for 
SU(N): 



d[Aj] 
d[kA x \ 

?4[fcAi] 



N(N - 1) • • • (N - j + 1) 
P- 

N(N + 1) • • • (N + k - 1) 
k\ 

(N - 4)! N- j 
N! (7-1)! 

(N-l)! N + k 
(N + 3)!(Jfc-l)! 



1 < / < N - 1 
k>l 



Q~[[(N-r + l)b 2 

k 

K k Y[(N + r-l)b 2 



d[Aj] 

d[kAi\ 

UAj] 

UkA,] 



l<j< 



N-3 



j\(N-j)\ ~ J - 2 

N + 2k-2(N + 2k-3)l 
k\ (N - 2)1 

(N - 4)! 



k> 1 



(;-l)!(N-;-l)! 

(N-2 + 2k)(N -2 + k)\ 
' (Jfc-l)!(N + 2)! X 
[N 2 - 3N + 8 + 6k(N -2 + k)\ b 2 



for N odd: 

d[A i v_ 1 ]=2^ 1 , J 4 [Aam] = -2 
for N even: 

d[A K ] = 2 t ¥, J 4 [An] = 



- N-10 

-2~ 



Sp(N): 



d[Aj] 
d[fcAi] 

r 4 [A ; ] 

l4[*Al] 



1 < 7 < N/2 



N + 2-2j (N + 1)! 

: yi (N + 2-7)! 
s CN±*zM k>1 

kl(N-l)l 

(N + 2- 2j)(N - 3)! 

: (i-iy.(N-j + iy. x 

X [N 2 + 3N + 8 - 6j(N + 2 - ;')] b 2 
(N + k)\ 



(jfc - 1)!(N + 3)! 



In particular, we list the coefficients for the relevant 
groups in this paper, where b was taken as 1: 



r\h[r] 


SU(N) 


SO(N) 


Sp(N) 


G 


2Nb 2 


(N - 8)b 2 


(N + 8)b 2 


E 


b 2 


b 2 


b 2 


m 


(N + 8)b 2 


(N + 8)b 2 


(N + 8)b 2 


B 


(N - 8)b 2 


(N - 8)b 2 


(N - 8)b 2 



1-1 
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